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Abstract 

We develop a matrix model to describe bilayered quantum Hall fluids for a series of 
filling factors. Considering two coupling layers, and starting from a corresponding action, 
we construct its vacuum configuration at v = q%K^ qj, where Kij is a 2 x 2 matrix and 
qi is a vector. Our model allows us to reproduce several well-known wave functions. We 
show that the wave function ^{m,m,n) constructed years ago by Yoshioka, MacDonald 
and Girvin for the fractional quantum Hall effect at filling factor m ^_ n and in particular 
\P(3,3,l) at filling I can be obtained from our vacuum configuration. The unpolarized 
Halperin wave function and especially that for the fractional quantum Hall state at filling 
factor | can also be recovered from our approach. Generalization to more than 2 layers 
is straightforward. 
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1 Introduction 



The quantum Hall (QH) effect has bred much beautiful theory. Indeed, Laughlin's wave func- 
tions [TJ are good wave functions for describing the fractional quantum Hall effect (FQHE) [2j at 
filling factor v — — , where m is an odd integer. For other filling factors several attempts have 
been suggested to extend Laughlin's theory by adopting different approaches and assumptions. 
In particular, Halperin p] proposed a family of generalized Laughlin wave functions that could 
incorporate reversed spins. In fact a candidate for an unpolarized wave function at filling factor 
| was given. Subsequently, Yoshioka, MacDonald and Girvin || generalized the Laughlin wave 
functions to those of the bilayered QH systems and derived that corresponding to the v = | 
state. Moreover, other theories have been elaborated and led to understand the observed values 
of v, in particular v = | [5J as well as others jHj- 

The first experimental indications for an unpolarized ground-state spin configuration in the 
FQHE came with the discovery of the v = | state [7j and later the v = | state 0. More 
compelling evidence for novel spin phenomena in the FQHE was subsequently reported j^j. 
On the other hand, it was shown experimentally that multi-layer systems also exhibit the 
FQHE [HI] ■ in fact, several filling factors have been observed, for instance the v — ~ state ^T] 
and i/= f,^,-" LIS- 

Recently, Susskind [T3] proposed a novel method to investigate the FQHE. He claimed that 
the non-commutative Chern-Simons theory (NCCS) at level k is exactly equivalent to Laughlin's 
theory at the filling factor us = \- He formulated his approach as a matrix theory similar to that 
describing DO-branes in string theory. However, Susskind's theory is an alternative approach 
to the FQHE which so far has not produced anything new but has just recovered the Laughlin 
approach by adopting a new formalism. Nevertheless it remains a new way of thinking and is 
worth studying in the hope that it will bring new results in the future. 

Although the proposed matrix model seems to reproduce the basic features of the Laughlin 
QH droplets, still some problems remain to be solved. Indeed, Susskind's approach is valid 
only for infinite matrices and also shows an anomaly for k = 0. To solve these problems, 
Polychronakos J3] introduced a boundary term to Susskind matrix model. He proposed a 
finite matrix model as a regularized version of the NCCS theory. It allowed him to find a 
quantum correction to v^, where k is shifted to k + 1 and the filling factor became z/ps = 
As another consequence, he pointed out that his matrix model is equivalent to the Calogero 
model [To] . 

Sometimes later, observing that the Laughlin wave functions can be mapped onto many- 
body wave functions of the harmonic oscillator, Hellerman and Van Raamsdonk ^0] built a 
complete minimal basis of wave functions of the theory at arbitrary level k and rank M, see 
also J7j. Other investigations about the relation between NCCS and Laughlin fluids can be 
found in [THJ EH]- Subsequently, the Susskind model and its regularized version introduced 
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by Polychronakos was extended to FQH states that are not of Laughlin type: a multicompo- 
nent Chern-Simons approach was introduced [201 and another proposal based on the Haldane 
hierarchy [2Tj was developed [22| 123] . 

Despite the progress in the study of the FQH fluids in the framework of NCCS matrix 
model, several open questions remain which have not been addressed so far. One of these 
questions concerns the wave functions that are not of Laughlin type. In fact there are many 
wave functions, that have been constructed years ago, e.g. by Yoshioka et ai, Halperin, • • -, 
but cannot be recovered by what is developed so far. 

In what follows we propose a matrix model to investigate the possibility to obtain two of 
those wave functions. This can be done by extending the Susskind-Polychronakos model to 
deal with the QH fluids at the filling factor [21] 



where is an iV x A" matrix and is a vector. The basic idea is to consider several Susskind- 
Polychronakos systems, let us say M systems, with an interaction between them and suppose 
that all systems possess the same number of particles. In the QHE language, this picture is 
equivalent to considering multi-layered systems. Without loss of generality, we fix M = 2, but 
as we will see later our analysis can directly be extended to the generic case M > 3. 

We start by writing down an appropriate action as a sum of two terms, for the free and 
the interacting part. Subsequently, we derive the corresponding Hamiltonian, which of course 
contains an interaction. Using a unitary transformation, we show that this Hamiltonian can be 
transformed to a diagonalized one. Next, we determine the vacuum configuration that allows 
us to recover two different states. Indeed, we show that how the Yoshioka-MacDonald-Girvin 
wave functions at the filling factor v = -^-^ can be obtained from our model and in particular 
that describing the FQHE at v — ^. Moreover, the unpolarized Halperin wave functions will 
be derived and especially that corresponding to v — | state. 

In section 2 we recall briefly the NCCS matrix model description of the Laughlin fluid. In 
section 3, we propose an action describing a system of two layers, we derive the Gauss law 
constraint as well as the equations of motion for the different variables. A quantum mechanical 
analysis will be the subject of next section, where we develop a Hamiltonian that corresponds 
to the system under consideration. Under rotation, we define a set of matrices of harmonic- 
oscillator operators to diagonalize the system. In section 5, we build the vacuum configuration 
that satisfies the constraint. A link with literature will be discussed in the last section where 
the two wave functions mentioned above will be recovered. We conclude our paper by putting 
some questions to be investigated in forthcoming works. 
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2 Chern-Simons matrix model 



Starting from the matrix formulation of a two-dimensional system with a large number of 
electrons in the presence of a perpendicular magnetic field B, Susskind ^Hl showed that the 
resulting effective theory is a non-commutative U(l) Chern-Simons gauge theory at level k = 
B6. As a consequence, he found a relation 

1 



2n6 



which links the non-commutivity parameter 9 to the density of electrons p. By using the 
definition of the filling factor 

y = ^ (3) 

in the system of units (H,e,c), it is easily seen that the fraction v can be written in terms of 
the parameter 6 as 

" = w (4) 

This beautiful relation is one of the interesting results obtained recently by Susskind in dealing 
with the FQH fluids. 

Moreover, by exploring the possibility to develop a consistent finite matrix model for the 
description of the FQH droplet, Polychronakos [H] suggested to include a new field into the 
Susskind model. The proposed action is given by 

S = j dt | Tr {e ab (± a + i[A , Xj) X b + 26 A, - uX^ + ^ - A ^j (5) 

where X a , a = 1,2 are N x N matrices and ip is a complex iV-vector, and e 12 = — e 21 = 1, 
e aa = 0. The action is invariant under the gauge group U(N) and the matrix model variables 
transform as 

X a ^UX a U-\ iP^Utp. (6) 

The equation of motion for A leads to the Gauss law constraint 

G = -iB [Xi, X 2 ] + V^ f - B9 = 0. (7) 
The trace of this equation gives 

ft if> = NB6. (8) 

Upon quantization the matrix elements of X a and the components of i/j become operators, 
obeying the commutation relations 



ipi,ftj 

[(Xi)jj, (X 2 )fe/] = jjSu 5jk- 



(9) 



3 



The Hamiltonian can be obtained from (jSJ) as 



f N 2 \ 



(10) 



where the N x N matrix of harmonic-oscillator operators is defined by 



A nm = ^(X 1 + iX 2 ) nm . (11) 

The corresponding wave function is ^H] 

\k) = [e^-^l i (^Ut) J2 ...(^t^-i )jjv ] /£ | ) (12) 

where the vacuum |0) is annihilated by A's and -0's and e is the fully antisymmetric tensor. 
This is a physical state and therefore satisfies the relation 

G\k) = 0. (13) 

It is similar to the Laughlin wave function 1J at the filling factor 

v = - . 14 

k+ 1 v ; 



Subsequently, one of us and others (221 EHj generalized the above results to any filling factor 
which can be expressed as 

ffata = T + t ( 15 ) 
k\ k 2 



and in particular to level two of the Haldane hierarchy 

P2 



(16) 



^ PlP2 ~ 1 

by setting 

h = Pi, k 2 = Pi(pip 2 - 1) (17) 
where pi is an odd and p 2 is an even integer. 

3 Two coupling matrices model 

We consider two systems with a total number of particles Mi + M 2 which interact with each 
other. Such systems can be seen like two coupling layers i containing Mj particles. The 
appropriate action to describe the FQH fluids of the whole system at filling factor ((TJ), is given 
by 

S = I ^Et^fl 1 ^^ 4 X^])X^ + 29A -U J (X^) 2 } 

J 3 

+V ,0-)t ^0) _ Aq ^ + f dt Kn j^i Tr (xWxW) + ^ (1 V (2) } (18) 
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which involves two copies of the single- layer action (J3J) forming the free part. It also contains 
an interacting part, where the scalar K\ 2 is playing the rule of a coupling parameter between 
the layers 1 and 2. The ratio ^f- is basically the magnetic field B. 

It is clear that for K\ 2 = 0, the total system becomes decoupling. Note that as far as 
the total action is concerned, the full gauge symmetry is U(M{) x U(M 2 ). The matrix model 
variables transform under this invariance as 

X« -> UX&U- 1 , -> U^. (19) 

Compared to the original matrix model, there is the potential term 

V = £ l^ Tr - (20) 

analogous to the potential of two coupled harmonic oscillators 5^1 in two-dimensional space. 
This provides a Hamiltonian for the theory. 

The Gauss law constraint can be obtained by evaluating the equation of motion for Aq. In 
our case it reads 



Q = -iK 



11 



x£\x®} - zK 22 [xfUf 1 + (^ (1 V (1)t + V; (2 V 2)t - tfn - K 22 ) = (21) 



where its trace gives 

^(i)t^(i) + ^(2)t^(2) = M x K n + M 2 K 22 . (22) 
Other equations of motion can also be calculated. For the X's we get 

Kue^X^ + KnunX® + K 12 u 12 xP = 
K 22 e ab X { a ] + K 22 u; 2 xi 2) + K 12 u l2 x£ ] = 

while for the if)'s we obtain 

^(Dt + Kl ^(2) = 

^(2)t + Kl2 ^ = 0. 1 J 

Of course the last set of equations shows a difference with respect to the decoupled case. It can 
be solved by using a unitary transformation. 



4 Hamiltonian formalism 

Let us now consider the proposed model quantum mechanically. We proceed by determining 
the total Hamiltonian, which describes the system under consideration. It can be obtained 
from the relation 

r)T 

U = X—r - L (25) 

dx v ; 
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where J~ defines the conjugate momentum. This leads to a Hamiltonian as the sum of the free 
and the interacting parts as 

H = E lf^ Tr - t 1 " 1271 w i)x « 2) ) (26) 

which is nothing but the confining potential (|2(J|) . This means that the kinetic energy is negli- 
gible compared to V. 

It is clear that this form of 7i cannot be diagonalized directly. Nevertheless, 7i can be 
transformed to another, factorizing Hamiltonian H! . Probably the best way to do this is to 
perform a rotation by a mixing angle ip, of the X's to new matrices 

y a (1) = xPcosf -Xi 2) sinf 

yi 2) = X«sinf + jd 2) cosf. ( ' 7j 

It can easily be checked that this rotation is a unitary transformation. Inserting (J27|) into (|26J1 . 
one can show that Tt transform to 

^ = «Tr(yW) 2 + /3Tr(ri 2 )) 2 (28) 

if the rotating angle satisfies the relation 

K 12 u 12 , 

tan tp = — - . (29) 

KnUi - K 22 uj 2 

The parameters a and (3 are given by 

a = | (KnCJi cos 2 | + sin 2 f - \Ki 2 u X2 sin </?) 

/? = I (i^n^i sin 2 I + K 22 u 2 cos 2 | + \K X2 uj 12 sin 99) . 

To diagonalize 7i\ we define two couples of creation and annihilation matrices of harmonic 
oscillator operators, 

2 )nm 
2 )nm- 



(1) 

nm — 




(2) 

nm — 


v^l (2) - 


;ions 






^nm 


r (l)t 

n'm' 




>(2) 




— 



(32) 

Oij'Oi'j, 

while all others commutators vanish. Now 7i' can be rewritten as 

W' = £ (2M X + Ml) + ^ (2M 2 + M 2 ) (33) 

_ _ 

where the number operators 

J v l\ — Z_m m=l ^ mn ^"nm 

m 2 = yJIU c ^ c \ 

are counting the M\ and M 2 particles. Thus under the unitary transformation the system 
became decoupling. 



M 2 n {2)] n {2) 

(34) 
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5 Ground-state wave functions 



To begin we emphasize a difference between the ground state of two coupled harmonic oscillators 
in terms of the coordinates Xi and that in terms of their mapped representations yi. The wave 
function 

t/j (y) ~ exp {-ayf - /3y%} (35) 

is separable in the variables y\ and y 2 . However, for the variables X\ and x 2 , the wave func- 
tion (f3~K|) reads 



ijj (x) 



exp 



-a [ x\ cos — — X2 sin 



(p\ 2 ( (p ( P\ 2 1 

in - J - (3 [xi sin - + x 2 cos - J > . 



(36) 



Next, we will see how these ground states can be extended to the matrix model formalism. 
We begin to determine that for the matrices Y. By transforming the Gauss law constraint to 
the variables Y, i.e. 



K u cos 2 - + K 22 sin 2 - 



{K n -K 22 )sm V { y} 1] ,Y^ 2) 



y(l) V W 

I I > 2 2 



+ ( K n sin 2 % + K 22 cos 2 % 



v m v (2) 

M > J 2 



+ 



lf\ F«] } = ifl (K n + K 22 - ^tyM - 0( 2 y 2 )t) (37) 
where the Polychronakos fields are also rotated to new fields 



0(1) = ^(1) CQS | _ ^(2) gin | 
0(2) = ^,(1) gin | + ^(2) cog |_ 

For simplicity let us fix Ku = K 22 = K, then (p?7|) becomes 



(3? 



1 1 I 2 2 



+ 



V( 2 ) V( 2 ) 
M 5 2 2 



2/7x6; I 1 - -L^tywt 
2iY 



2iY 



0( 2 )0( 2 )t 



(39) 



Now it is clear that the ground state is simply a tensor product between those states corre- 
sponding to each layer 



\K) 



^••^4 i 1 )t(0t(i) C '«t) i2 . . . (^W^W-i). 

cii-« 2 0(. 2 )t(0(2)tC(2)t) . 2 . . . (0(2)t C (2) t M 2 -l 



Ml 



A' 



A 



|0>. 



(40) 



The ground state (J4(J|) can be mapped in terms of the operators of the matrices X by 
expressing the matrices C of harmonic-oscillator operators in terms of those corresponding to 
the matrices X. Using (|27jl one can show that (j31j) takes the form 



a 



(2) 
nm 



f (A« cosf - sin|) 7 
7§(A« sin I + A( 2 ) cosf) ] 



(41) 
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where the operators 



A 0) 

■rir, 



fv(l) ,:y( 



commute: 



4(1) .Wt 

4(2) ,(2)f 



SifSi'j. 



(42) 



(43) 



Inserting (j3E|) and (}4~Tj) in (j4*U|) . we obtain 

\k)-- 



2 1 / ii 









»Afl " 



h—3M 2 ^Wtsin^ + ^tcos^) ••• 
\ 2 2 / ji 

I ty,(Dt sin £ + ^(2)t cos 1) (^(Dt sin £ + A (2)t cos § ) M2_1 



|0>. 



(44) 



In what follows, we proceed without the use of the unitary transformation to construct the 
wave function |$) describing the system of M\ + M 2 electrons at filling factor (0). One has to 
realize a physical state |$) that satisfies the Gauss law constraint (|2~Tj) 



g\$) = 



(45) 



and allows us to establish a link with two well-known wave functions. May be the best way to 
do this is to define two operators 

A = A^®A^ 

4 = v>w ® ^ (2) (46) 

where ® is the tensor product. Using these matrices of harmonic-oscillator operators, we build 
a vacuum configuration 

Ku-K l2 

1>M\ 



(^ (1)t ^ (1)t ) 4 • • • (^ (1)tM1 - 1 ) ! 



ix~4u^®\ (^A^) n ■ ■ ■ (^(2)t A (2)tM 2 -l) 



which satisfies the Gauss law constraint (j21j) and therefore we have 

(^ (1 V (1)t + ^ (2 V (2)t - M x K n - M 2 K 22 ) |$) = 0. 
Novel about this vacuum configuration is that one can interpret the term 



K22—K12 

K 12 



|0>. 



,fci...fcM-,+M 9 



(47) 



(4f 



(49) 
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as an inter-layer correlation. In conclusion, our configuration could well be a good ansatz for 
the ground states of double-layered FQH fluids in the formalism of the NCCS matrix model. 
This will be clarified in the next section. 



6 Link with literature 



Here we show how the Yoshioka-MacDonald-Girvin and Halperin wave functions describing, re- 
spectively, the double-layer and the unpolarized QH systems can be recovered from our vacuum 
configuration (JJ7J). 

Before starting, we note that for any iV-dimensional vector ij)< and N x N matrix A\ the 
expression of the form 

Fty\ At) = ^-^(D ^(i)t A (i)t) ^ . . . (^(Dt^DtA-i^ 

has a one-to-one correspondence to the polynomial 

f{z) = £ n.., Nz O i ... z N-l_ 

Now our task can be done by defining a new complex variable 

'% 

,( 2 ) 



in 



: , for i = 1, N 



Z i-N 



(50) 



(51) 



(52) 



fori = iV + l,...,2iV 

assuming that the particle numbers are equal, Mi = M 2 = N, and recalling the Vandermonde 
determinant: 

(53) 



n 

Kj 



Zi — Zj) = det ( z- 



In terms of the complex coordinates, reads 



n , (1) xAT-l 
Z i N 



N-l 



Ku-K 12 



I<22 — K 1 2 



It can be written in standard form as 



(i) _ ^y 11 

i<j i<j 



k 1 ...k 2N /-0 . . . X2JV-1 

n 



1^12 



*i ) 



and now the inter-layer correlation is 



n 



i,0 



(54) 



(55) 



(56) 



Next, we will give two different applications of 1)55)1 . 
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6.1 YMG wave functions 

Considering the two layers and treating them as additional degrees of freedom, the v = | 
state was predicted by Yoshioka, MacDonald and Girvin [I]. They made a straightforward 
generalization of the Laughlin wave functions to those with the filling factor 

(57) 



m + n 

where m and n are integers. This can be obtained from our analysis by taking 

K=i m " ), 9 =(l -l) (58) 

leading to the wave function 

-n(«p- -f ) "' n if - -f *) " n (*."'- -f ') " *»• m 

i<j i<j i,j 

Choosing m = 3 and n — 1, we recover the FQHE = | state corresponding to 

= n - #') 3 n (*r - f o'n - »- < 6 °) 

«<j i<j i,j 

6.2 Halperin wave functions 

Another interesting result can be obtained. In the Halperin picture in the context of single- 
layered unpolarized QH systems, the labels 1 and 2 can be considered as an analogue of spin. 
Following this idea, our bilayered system can be seen as mixing layers of particles with spin up 
and spin down. 

As a consequence, we obtain for m = 3 and n = 2 the unpolarized Halperin wave function 
with the filling factor | as 

*w* = n (*. (,) - #') 3 n o. ,2 > - # o 3 n (« w - *- < 6i > 

i<j i<j i,j 

This can be seen as a wave function of a system of N particles with spin parallel and another 
N particles with spin antiparallel to the external magnetic field. 

7 Conclusion 

We have developed a matrix model to describe bilayered QH systems at the filling factor 
v = qiK^qj. The basic idea was to use two coupled harmonic-oscillators in a similar fashion 
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as done by Susskind and Polychronakos. Our model is a generalization of their model and of 
course reproduces its basic features by taking the coupling parameter K 12 to be zero. 

Starting from an appropriate action we derived the equations of motion for the different 
matrix model variables. The corresponding Hamiltonian was obtained as the sum of free and 
interacting terms. A unitary transformation, more precisely a rotation around an angle ip, led 
to a factorizing Hamiltonian. 

Next, we have constructed the ground states of the system in two different ways. The 
first was based on the unitary transformation and from the ground state after rotation we 
have derived that before rotating the system. The second was performed directly in terms of 
a combination of the matrices of harmonic-oscillator operators of two layers. The obtained 
vacuum configuration involved three different quantities where one describes the inter-layer 
interaction. 

Subsequently, we have investigated the link between our second wave function and two 
others from literature. After projecting the vacuum configuration on the complex plane and 
using the Vandermonde determinant, we have shown how the Yoshioka-MacDonald-Girvin 
wave function with the filling factor v = can be obtained from our model, in particular 
that corresponding to the v — | state. Likewise, we have recovered the unpolarized Halperin 
wave function, especially that for the v = | state. 

The case we have studied is in fact just a particular case of more general FQH states where 
the fluid droplet is assumed to consist of several coupled branches, say M branches. M = 1 is 
the Laughlin (Susskind-Polychronakos) model, M = 2 is the model we have discussed here and 
M > 3 is the generic case, which can be seen as a straightforward generalization of our case. 

Of course still some important questions remain to be answered e.g. about the fractional 
charge and statistics of the particles and how to describe them in terms of the proposed model. 
Another interesting question is related to the link between our model and Calogero and super- 
Calogero models. We will return to these issues and related matter in future. 

We close this section by noting that our model will be investigated in the forthcoming 
work [2H1 for the case of a single layer. Basically, we will consider the Laughlin liquids in a 
confining potential that is not of parabolic type and see how this affects the basic features of 
the Susskind-Polychronakos model. 
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